Vehicle loading based vehicle dynamic and safety related characteristic adjusting system

ABSTRACT

A method of controlling a controllable chassis system or a safety system ( 44 ) for a vehicle ( 10 ) includes determining an added mass placed on the vehicle and relative to a known vehicle mass. A vehicle characteristic is adjusted in response to the added mass. A control system ( 18 ) for an automotive vehicle ( 10 ) includes a sensor ( 20, 28 - 42 ) that generates a signal. A controller ( 26 ) determines added mass on the vehicle ( 10 ) in response to the signal and adjusts a vehicle characteristic in response to the added mass.

RELATED APPLICATIONS

This application is a divisional of Application Ser. No. 11/536,105, filed on Sep. 28, 2006, now U.S. Pat. No. 7,826,948, incorporated by reference herein. U.S. patent application Ser. No. 11/536,105 is a continuation-in-part of U.S. patent application Ser. No. 10/966,395, filed on Oct. 15, 2004 and now U.S. Pat. No. 7,668,645 also incorporated by reference herein.

TECHNICAL FIELD

The present invention relates generally to vehicle control systems and vehicle motion sensing systems. More specifically, it relates to techniques for adjusting ride comfort, handling and safety related characteristics of a vehicle in response to vehicle loading variation and loading distribution variations.

BACKGROUND

Vehicle roll stability control (RSC) schemes are used to address the issue of friction-induced rollovers. RSC systems include a variety of sensors that sense vehicle states and a controller that controls a distributed brake pressure to reduce a tire lateral force such that the net moment of the vehicle is counter to the roll direction.

During an event causing the vehicle to roll, the vehicle body is subject to a roll moment due to the coupling of the lateral tire force and the lateral acceleration applied to the center of gravity of the vehicle body. This roll moment causes suspension height variation, which in turn results in a vehicle relative roll angle (also called chassis roll angle or suspension roll angle). The relative roll angle is an important variable that is used as an input to the activation criteria of RSC and to construct the feedback brake pressure control command. The sum of such a chassis roll angle and the roll angle between wheel axle and the road surface (called wheel departure angle) provides the roll angle between the vehicle body and the average road surface, which is another important variable feeding back to the roll stability control module.

Such a chassis roll angle can be calculated as in U.S. Pat. No. 6,556,908 using the lateral acceleration of the center of gravity of the vehicle body and the roll angular acceleration, together with vehicle parameters such as the sprung mass, the vehicle body roll moment of inertia, the roll stiffness and damping ratio of the suspensions and the anti-roll-bars, and the distance between the center of gravity of the vehicle body and the floor of the vehicle body. The disclosure of U.S. Pat. No. 6,556,908 is hereby incorporated by reference.

One problem with using these parameters in the computation is that the aforementioned relative roll angle may vary with the vehicle operating conditions. For example, a 150 pound roof loading for a typical SUV with a curb weight of 5000 pounds may cause more than 30% error in relative roll angle calculations if computed assuming no roof load. From the vehicle mass point of view, although a 150 pound roof loading accounts for only a 3% mass variation over the vehicle curb weight, it could account for a 30% error in the chassis roll computation, which is ten times larger. If the above parameters are fixed at certain nominal values in the RSC system, it is conceivable that optimal control performance may not be achieved under a different loading condition. For example, if the relative roll angle is computed with nominal vehicle loading condition assumptions, without considering roof loading, the relative roll angle may be under estimated for vehicles with roof loadings, which results in a reduced control. That is, the control system may not be as effective as desired. On the other hand, if the relative roll angle is computed with maximum roof loading, it may be over estimated for vehicles without roof loadings causing unintended control. That is, the control system may become too sensitive or intrusive. Therefore, in order to improve the overall performance of the RSC system, it may be desirable to estimate and update the vehicle parameters periodically or adaptively in real time based on the detected roof loading.

Certain schemes for obtaining vehicle parameters have been disclosed. For example, in U.S. Pat. No. 4,548,079, a method is disclosed for determining vehicle mass directly using engine torque and vehicle acceleration. Similarly, in U.S. Pat. No. 5,490,063, push force is determined from the driveline torque and gear ratio to obtain vehicle mass. In U.S. Pat. No. 6,167,357, instead of calculating vehicle mass directly, a recursive least square (RLS) algorithm is proposed to estimate both vehicle mass and aerodynamic coefficient online. The latter method is considered to be more reliable since it recursively adjusts for estimation error of the previous estimates. Furthermore, the use of vehicle acceleration, which is usually very noisy, is avoided. The mass estimation schemes proposed in the above-cited patents may not accurately indicate changes to parameters that impact the roll dynamics of the vehicle. For example, a 150 pound roof loading on a 5000 pound SUV, i.e., 3% mass change, might be undetectable in the above schemes due to the potential error in the engine torque, which usually is much larger than 3%. Other error sources include the tire rolling radius change due to tire pressure drop and due to the vehicle loading variations, the vehicle drag and the offset or uncertainty in the longitudinal accelerometer.

The above schemes focus mainly on large mass variations, which may have significant influences on the vehicle longitudinal dynamics and vehicle fuel consumption. They do not differentiate when the vehicle mass change is due to a floor loading or due to a roof loading. However, the roof loading causes much more significant roll motion parameter changes than does the same amount of floor loading. That is, there is a need to detect not only the amount of loading (maybe much smaller than that can be detected by the existing method), but also the location of the loading (the vertical and longitudinal distance of the loading with respect to the vehicle floor or the center of gravity of the vehicle body, for example).

Furthermore, the other parameters that affect vehicle body roll and lateral dynamics, such as the roll stiffness and damping in the suspensions and the total center of gravity height of the vehicle body with respect to the vehicle floor, the roll moment of inertia, should also be taken into account.

With the advance in electronic controls, currently, some actuators designed for vehicle's ride and handling performances are now finding their ways to be utilized for stability control and safety. For example, suspension controls might be used for achieving certain yaw stability control function through individual suspension controls. The same is true for the actuators designed for stability control are utilized for vehicle's comfort and ride handling performances. For example, the selective braking and engine torque reduction used in stability control might be utilized for adjusting the vehicle's response during non-stability event such that the driver is experiencing a smooth and less intrusive driving. As another example, the control suspensions to achieve the vehicle's ride comfort and handling performance usually adjust only the suspension parameters to counteract these road disturbances such as gravel and potholes and they are not necessarily designed to account for vehicle loading, which, as stated above, can affect vehicle response too.

Therefore, there is a need for a technique that can detect or adaptively update the vehicle parameters, including vehicle loading and loading location, the roll moment of inertia, in order to refine and improve a vehicle control system relating to both the stability control functions including roll stability control function during vehicle stability control events and the ride comfort and handling functions during normal driving events.

SUMMARY OF THE INVENTION

It is one object of an embodiment of the invention to provide a system for a quantitative characterization of the chassis roll that may be used in conjunction with the stability control system of a vehicle to determine accurate vehicle roll trending for vehicles with varying inertia parameters due to loading conditions, and adaptively adjust the computation of the chassis roll angle and/or adaptively adjust the activation criteria for the stability control function. More specifically, the present invention seeks to determine the roll condition in response to an additional mass and/or height of the mass, that is, to determine if a small mass change may significantly impact the roll dynamics of the vehicle, as is the case of a roof load. The control system then can make a determination with respect to adjusting vehicle control command. In one embodiment, the sensitivity and profile of the control commands of various actuators and/or the response thereof are adjusted.

In another embodiment, a method of controlling a vehicle dynamics includes determining an added mass placed on the vehicle and relative to a known vehicle mass. A vehicle characteristic, such as the ride, comfort and handling, is adjusted in response to the added mass.

In a further embodiment, a control system for an automotive vehicle includes a sensor that generates a signal. A controller determines the added mass on the vehicle in response to the signal and adjusts at least one vehicle characteristic in response to the added mass.

In yet another embodiment of the invention, a method of controlling a vehicle is provided and includes determining a roll gradient, a roll acceleration coefficient, and/or a roll rate coefficient. That is, the roll acceleration coefficient, roll gradient and the roll rate coefficient may be used simultaneously or in various combinations to control the vehicle dynamics.

In still another embodiment of the invention, a method of controlling a control device for a vehicle is provided that includes determining a composite parameter called roll gradient, determining another composite parameter called roll acceleration coefficient, determining an added mass and a height of the added mass from the roll gradient and the roll acceleration coefficient, and controlling the control system in response to the added mass and height.

In another embodiment of the invention, a control system for an automotive vehicle is provided and includes a roll, rate sensor generating a roll rate signal, a lateral acceleration sensor generating a lateral acceleration signal, a longitudinal acceleration sensor generating a longitudinal acceleration signal, and a yaw rate sensor generating a yaw rate signal. The control system and the sensors are coupled to a controller. The controller determines an added mass and a height of the added mass from the roll rate, the lateral acceleration, the longitudinal acceleration and the yaw rate and controlling the control system in response to the added mass and height.

In another embodiment of the invention, a method of controlling a control device for a vehicle is provided and includes determining different roll conditions differentiated through determining an added mass and a height of the added mass, and controlling the control system in response to the roll condition due to an aggressive driving maneuver and enhanced by the added mass at the given height.

One advantage of the invention is that various vehicle loading conditions can be detected and determined very accurately. Such loading conditions (for example, roof loadings) are usually a fraction of the total vehicle mass (for example, 3% to 6% of the vehicle curb weight), are difficult to be detected by using the methods involving longitudinal vehicle dynamics, and may cause large (for example, more than 30%) changes in the chassis roll angle.

Another advantage of the present invention is that it may be utilized alone or in conjunction with or as part of various vehicle control systems to improve vehicle ride, handling and safety performances for various normal and abnormal driving conditions.

Other advantages and features of the present invention will become apparent when viewed in light of the detailed description when taken in conjunction with the attached drawings and appended claims.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagrammatic view of a vehicle with variable vectors and coordinate frames according to the present invention.

FIG. 2 is a block diagram of a stability system according to the present invention.

FIG. 3 is a front view of an automotive vehicle illustrating various angles according to the present invention.

FIG. 4 is a side view of an automotive vehicle illustrating various variables thereon.

FIG. 5 is a top view of an automotive vehicle having variables used in the following calculations thereon.

FIG. 6 is a block diagrammatic view of a controller according to a first embodiment of a least square method of the present invention.

FIG. 6A is a block diagrammatic view of a controller according to a second embodiment of a least square method of the present invention.

FIG. 7 is a logic flow diagram illustrating a method of operating a safety system in accordance with multiple embodiments of the present invention.

DETAILED DESCRIPTION

In the following figures, the same reference numerals will be used to identify the same components. The present invention may be used in conjunction with a roll stability control system for a vehicle. The present invention may also be used with the other chassis control system such as a controllable suspension system, a controllable steering system, an active roll bar. The present invention will be discussed below in terms of several embodiments relating to an automotive vehicle moving in a three-dimensional road terrain. Some embodiments of the present invention are described with respect to determining an added mass and height of the mass. As will be described below the added mass and height may not be directly determined, rather by adaptively updating a roll gradient value, roll acceleration coefficient, roll moment of inertia and/or a roll rate coefficient, the effects of added mass and the height may be included in those values. Such values may also be referred to as an “adaptive” roll gradient, an “adaptive” roll acceleration coefficient, an “adaptive” roll moment of inertia and an “adaptive” roll rate coefficient. The various roll parameters may be referred to as adaptively determined, meaning that such values may change due to a changing mass or load over time. That is, such values are not fixed as has been done in previous systems.

It should also be noted that while one of the roll gradient, roll acceleration coefficient and roll rate coefficient may be adaptively determined, the other two values may not be adaptive.

Referring to FIG. 1, an automotive vehicle 10 with a safety system of the present invention is illustrated with the various forces and moments thereon during a rollover condition. Vehicle 10 has front right (FR) and front left (FL) wheel/tires 12A and 12B and rear right (RR) wheel/tires 13A and rear left (RL) wheel/tires 13B, respectively. The vehicle 10 may also have a number of different types of front steering systems 14 a and rear steering systems 14 b, including having each of the front and rear wheels configured with a respective controllable actuator, the front and rear wheels having a conventional type system in which both of the front wheels are controlled together and both of the rear wheels are controlled together, a system having conventional front steering and independently controllable rear steering for each of the wheels, or vice versa. The front and rear wheel steering systems 14 a and 14 b may be or be part of an electronic power steering (EPS) or an active steering system. Generally, the vehicle has a weight represented as Mg at the center of gravity of the vehicle, where g=9.8 m/s² and M is the total mass of the vehicle.

As mentioned above, the system may also be used with chassis control systems and safety systems including active/semi-active suspension systems, anti-roll bar, or airbags or other safety devices deployed or activated upon sensing predetermined dynamic and safety conditions of the vehicle.

The sensing system 16 is coupled to a control system 18. The sensing system 16 may comprise many different sensors including the sensor set typically found in a roll stability control including lateral accelerometer, yaw rate sensor, steering angle sensor and wheel speed sensor which are equipped for a traditional yaw stability control system together with a roll rate sensor and a longitudinal accelerometer. The various sensors will be further described below. The sensors may also be used by the control system in various determinations such as to determine a lifting event, determine a height and position of a mass, etc. The wheel speed sensors 20 are mounted at each corner of the vehicle and generate signals corresponding to the rotational speed of each wheel. The rest of the sensors of sensing system 16 may be mounted directly on the center of gravity of the vehicle body, along the directions x,y and z shown in FIG. 1. As those skilled in the art will recognize, the frame from b₁,b₂ and b₃ is called a body frame 22, whose origin is located at the center of gravity of the car body, with the b₁ corresponding to the x axis pointing forward, b₂ corresponding to the y axis pointing off the driving side (to the left), and the b₃ corresponding to the z axis pointing upward. The angular rates of the car body are denoted about their respective axes as ω_(x) for the roll rate, ω_(y) for the pitch rate and ω_(z) for the yaw rate. Calculations may take place in an inertial frame 24 that may be derived from the body frame 22 as described below.

The angular rate sensors and the accelerometers may be mounted on the vehicle car body along the body frame directions b₁,b₂ and b₃ which are the x-y-z axes of the sprung mass of the vehicle.

The longitudinal acceleration sensor is mounted on the car body located at the center of gravity, with its sensing direction along b₁ axis, whose output is denoted as a_(x). The lateral acceleration sensor is mounted on the car body located at the center of gravity, with its sensing direction along b₂ axis, whose output is denoted as a_(y).

The other frame used in the following discussion includes the road frame, as depicted in FIG. 1. The road frame system r₁r₂r₃ is fixed on the driven road surface (which is yawing and translating with the vehicle), where the r₃ axis is along the average road normal direction computed from the normal directions of the four-tire/road contact patches.

In the following discussion, the Euler angles of the body frame b₁b₂b₃ with respect to the road frame r₁r₂r₃ are denoted as θ_(xbr) and θ_(ybr), which are also called the relative Euler angles (i.e., relative roll and relative pitch angles, respectively).

Referring now to FIG. 2, roll stability control system 18 is illustrated in further detail having a controller 26 used for receiving information from a number of sensors which may include a yaw rate sensor 28, a speed sensor 20, a lateral acceleration sensor 32, a vertical accelerometer sensor 33, a roll angular rate sensor 34, a steering wheel (hand wheel) angle sensor 35, a longitudinal acceleration sensor 36, a pitch rate sensor 37, steering angle (of the wheels or actuator) position sensor 38, suspension load sensor 40 and suspension position sensor 42. It should be noted that various combinations and sub-combinations of the sensors may be used.

Controller 26 may include a signal multiplexer 50 that is used to receive the signals from the sensors 28-42. The signal multiplexer 50 provides the signals to a wheel lift detector 52, a vehicle roll angle calculator 54, and to a roll stability control (RSC) feedback control command 56. Also, wheel lift detector 52 may be coupled to the vehicle roll angle calculator 54. The vehicle roll angle calculator 54 may also be coupled to the RSC feedback command 56. The RSC feedback command 56 may include a torque controller 57. Vehicle roll angle calculator 54 is described in U.S. Provisional Applications 60/400,376 and 60/400,172, and in U.S. patent application Ser. No. 10/459,697, the disclosures of which are incorporated herein by reference.

A loading detector 58 may also be included in controller 26. The loading detector 58 may be used to determine an additional mass of the vehicle and a distance of the mass as will be described below.

In one embodiment, the sensors are located at the center of gravity of the vehicle. Those skilled in the art will recognize that the sensor may also be located off the center of gravity and translated equivalently thereto.

Safety device 44 may control the position of the front right wheel actuator 46A, the front left wheel actuator 46B, the rear left wheel actuator 46C, and the right rear wheel actuator 46D. Although as described above, two or more of the actuators may be simultaneously controlled. For example, in a rack-and-pinion system, the two wheels coupled thereto are simultaneously controlled. Based on the inputs from sensors 28 through 42, controller 26 determines a roll condition and/or wheel lift and controls the steering position and/or braking of the wheels.

Safety device 44 may be coupled to a brake controller 60. Brake controller 60 controls the amount of brake torque at a front right brake 62 a, front left brake 62 b, rear left brake 62 c and a rear right brake 62 d. Other safety systems such as an antilock brake system 64, a yaw stability control System 66 and a traction control system 68 may also benefit from the knowledge of the additional mass and position of the mass detected by the loading detector 58. This information may be used to modify the control strategy used in computing the controlled brake force.

Speed sensor 20 may be one of a variety of speed sensors known to those skilled in the art. For example, a suitable speed sensor may include a sensor at every wheel that is averaged by controller 26. The controller may translate the wheel speeds into the speed of the vehicle. Yaw rate, steering angle, wheel speed and possibly a slip angle estimate at each wheel may be translated back to the speed of the vehicle at the center of gravity. Various other algorithms are known to those skilled in the art. Speed may also be obtained from a transmission sensor. For example, if speed is determined while speeding up or braking around a corner, the lowest or highest wheel speed may not be used because of its error. Also, a transmission sensor may be used to determine vehicle speed.

Load sensor 40 may be a load cell coupled to one or more suspension components. By measuring the stress, strain or weight on the load sensor a shifting of the load can be determined.

The center of gravity C is illustrated with nominal mass M. A roll axis is also illustrated at a distance D from the center of gravity. a_(y) is the lateral acceleration.

Referring now to FIG. 3, the relationship of the various angles of the vehicle 10 relative to the road surface 11 is illustrated. In the following, a reference road bank angle θ_(bank) is shown relative to the vehicle 10 on a road surface. The vehicle has a vehicle body 10 a and wheel axle 10 b. The wheel departure angle θ_(wdo) is the angle between the wheel axle and the road. The relative roll angle θ_(xr) is the angle between the wheel axle 10 b and the body 10 a. The global roll angle θ_(x) is the angle between the horizontal plane (e.g., at sea level) and the vehicle body 10 a.

Another angle of importance is the linear bank angle. The linear bank angle is a bank angle that is calculated more frequently (perhaps in every loop) by subtracting the relative roll angle generated from a linear roll dynamics of a vehicle (see U.S. Pat. No. 6,556,908 which is incorporated by reference herein), from the calculated global roll angle (as in U.S. Pat. No. 6,631,317 which is incorporated by reference herein). If all things were slowly changing without drifts, errors or the like, the linear bank angle and reference road bank angle terms would be equivalent.

Referring now to FIGS. 4 and 5, an automotive vehicle 10 is illustrated with various parameters illustrated thereon. A change in mass ΔM is illustrated relative to nominal center of gravity C₀. The center of gravity moves to C and the added mass ΔM is added thereto. The change in mass or load ΔM is positioned at a distance H above the load floor 80. The nominal center of gravity C₀ is located at position h above the load floor 80. The distance between the new center of gravity and the nominal center of gravity C₀ is ΔH.

The longitudinal acceleration is denoted by a_(x) whereas the longitudinal velocity is denoted ν_(x). The lateral acceleration and lateral velocity is denoted by a_(y) and ν_(y) respectively. The steering wheel angle is denoted by δ_(w). The wheelbase of the vehicle is denoted by the symbol WB.

Referring now to FIG. 6, controller 26 is illustrated in further detail. Controller 26 receives various sensor signals, for example, the roll rate, the lateral acceleration, the longitudinal acceleration and the yaw rate from the respective sensors. From the sensor signals, the loading height and amount of loading, the roll gradient, roll rate parameter, roll acceleration coefficient and roll inertia may be determined. These terms may be used to generate a vehicle roll angle compensation term and a roll stability control gain/threshold compensation term. Instead of determining the load and height of the load, an adaptive roll acceleration coefficient and an adaptive roll gradient that inherently have the load and height of the load therein may be determined. As mentioned above, these values in prior systems were fixed without regard to loading. Here they are adaptively determined to correspond to the loading condition. For example, these terms may increase the sensitivity of the system due to a higher center of mass. As will be further described below, the controller uses the roll rate signal and performs a derivative in box 84 to generate the roll acceleration signal {dot over (ω)}_(x), which is then filtered by the first filter 86 to generate the variable X, which is provided to a conditionality averaging box 88. The lateral acceleration signal is filtered in the second filter 90 and provided to the conditionality averaging box 88 illustrated as Y. The longitudinal acceleration is filtered in the third filter number 92 and provided to conditionality averaging box 88. As will be further described below, the conditional averaging box generates signals that are coupled to a least squares fit determination 98 which provides variables α and β. The α and β variables are provided to a loading/loading height detection block 100. The conditional averaging block and least squares fit block 98 illustrate a first conditional least square method illustrated by box 99. The loading/loading height detection block 100 generates a vehicle roll angle compensation term 102 and an RSC control gain/threshold compensation term 104.

Referring now to FIG. 7, the operation of the controller 26 is described in further detail. In step 102, the various sensors such as the roll rate sensor, the lateral acceleration sensor, longitudinal acceleration sensor and yaw rate sensor are read.

The roll signal for control is calculated as the (θ_(x)−θ_(refbank)), i.e., the subtraction of the reference bank angle from the global roll angle. The roll signal for control may be used to control the roll stability control system or other system or device.

As mentioned above, the various errors not limited to sensor signal drift (for example, due to environmental temperature variations), zero offset, sensor signal scaling factor (sensor signal sensitivity) nonlinearity, and numerical errors in computation, may enter into the various computations at various times. Thus, in certain situations the wheel departure angle or the reference bank angle may not be accurate. The following description describes how these values can be updated in response to wheel lift/wheel grounded values.

In step 104, the roll gradient, roll acceleration coefficient computation and a roll rate parameter are determined. In the following discussion, the roll gradient, the roll acceleration coefficient and the roll rate parameter are defined and then computed based on the measured signals and the computed signals.

There are two roll angle computations in the roll stability control: one is the relative roll angle as computed in U.S. Pat. No. 6,556,908 and the other is the global roll angle computation as computed in U.S. Pat. No. 6,631,317. U.S. Pat. Nos. 6,556,9087 and 6,631,317 are incorporated by reference herein. The relative roll angle is a function of the lateral acceleration a_(y), the roll acceleration {dot over (ω)}_(x), the roll rate ω_(x), the roll gradient α and the roll acceleration coefficient β and a roll rate coefficient η. That is, θ_(xr)=ƒ(a _(y),{dot over (ω)}_(x);α,β,η,ω_(x))  (1)

In z-transformation, θ_(xr) =αT _(roll)(z)a _(y) −βT _(roll)(z){dot over (ω)}_(x)  (2) where T_(roll)(z) is a filter which reflects the linear roll model of the vehicle, i.e., it includes the roll stiffness and damping effect due to suspension dynamics. Based on the roll dynamics, the following is true:

${\theta_{xr} + {\frac{D_{roll}}{K_{roll}}{\overset{.}{\theta}}_{xr}}} = {{\alpha\; a_{y}} + {\beta\;{\overset{.}{\omega}}_{x}}}$

Performing a Laplace transform to convert this equation into the frequency domain equation, leads to the following computation θ_(xr)(s)=T _(Roll)(s)(αa _(y)(x)+β{dot over (ω)}_(x)(s)) where the transfer function T_(roll) is

${T_{Roll}(s)} = \frac{1}{1 + {\eta\; s}}$ with the roll rate coefficient η is defined as

$\eta = \frac{D_{roll}}{K_{roll}}$

Using the algorithm shown in U.S. Pat. No. 6,631,317, the global roll angle θ_(x) (the roll angle of the vehicle body with respect to sea level) can be obtained through an anti-drift integration of the Euler roll angle velocity (which is different from the roll angular rate) and a roll reference signal through a steady-state recovery filter. On level ground when the vehicle does not have wheels in the air, such a calculated global roll angle is the same as the chassis roll angle when the vehicle roll dynamics is within its linear range (for example, the cornering acceleration of the vehicle is under 0.35 g).

Therefore, on level ground when the vehicle does not have wheels in the air and the vehicle has moderate cornering acceleration, the algorithm used to compute global roll angle as in U.S. Pat. No. 6,631,317 is a function of the roll rate ω_(x), the yaw rate ω_(z), the relative roll angle θ_(xr) and the relative pitch angle θ_(yr) θ_(x) =g(ω_(y),ω_(z),θ_(xr),θ_(yr))  (3)

In z-transformation, such a computation can be expressed as in the following format θ_(x) =T _(adi)(z)[ω_(x)+ω_(z)θ_(yr) ]+T _(ssc)(z)θ_(xr)  (4) where T_(adi)(z) is the so-called steady-state recovery filter, which is used to recover the good low frequency portion of the roll angle removed by the anti-drift-integration filter. Here the relative roll angle (chassis roll angle) θ_(xr) is used to recover the true steady state value of the roll angle.

Notice that on level ground if the vehicle is driven without wheel lifting and with moderate cornering acceleration, the relative roll angle and the global roll angle are equal. Therefore, from equations (1) and (3) or (2) and (4), the relationships between the roll gradient α, the roll acceleration coefficient β, and the roll rate coefficient η may be found.

Although the computation of the roll gradient α and the roll acceleration coefficient β is the main focus of this disclosure, the roll rate coefficient η can be similarly determined. It can be adaptively determined independent of the determination α and β, or simultaneously determined together with α and β. Likewise α and β can be independently determined.

Substituting (2) into (4), the following is obtained θ_(x) =T _(adi)(z)[ω_(x)+ω_(z)θ_(yr) ]+αT _(ssc)(z)T _(roll)(z)a _(y) −βT _(ssc)(z)T _(roll)(z){dot over (ω)}_(x)  (5) where the relative pitch θ_(yr) is a function of the longitudinal acceleration and the pitch acceleration as in the following (see U.S. Pat. No. 6,556,908 for detail) θ_(yr) =m(a _(x),{dot over (ω)}_(y);κ,λ)  (6) and which can be further written as θ_(yr) =κT _(pitch)(z)a _(x) +λT _(pitch)(z){dot over (ω)}_(y)  (7) where T_(pitch)(z) is a filter which reflects the linear pitch model of the vehicle, i.e., it includes the pitch stiffness and damping effect due to the vehicle suspension dynamics. T_(pitch)(z) is different for vehicle deceleration and acceleration, i.e.,

$\begin{matrix} {{T_{pitch}(z)} = \left\{ \begin{matrix} {T_{pitchacc}(z)} & {{{if}\mspace{14mu} a_{x}} > 0} \\ {T_{pitchdec}(z)} & {{{if}\mspace{14mu} a_{x}} \leq 0} \end{matrix} \right.} & (8) \end{matrix}$ and the same is true for the pitch gradient κ and the pitch acceleration coefficient λ

$\begin{matrix} {\kappa = \left\{ \begin{matrix} \kappa_{acc} & {{{if}\mspace{14mu} a_{x}} > 0} \\ \kappa_{dec} & {{{if}\mspace{14mu} a_{x}} \leq 0} \end{matrix} \right.} & (9) \\ {\lambda = \left\{ \begin{matrix} \lambda_{acc} & {{{if}\mspace{14mu} a_{x}} > 0} \\ \lambda_{dec} & {{{if}\mspace{14mu} a_{x}} \leq 0} \end{matrix} \right.} & \; \end{matrix}$

If the pitch angular rate ω_(y) is not available, the above relative pitch θ_(yr) can be approximated as in the following θ_(yr) ≈κT _(pitch)(z)a _(x)  (10) which is acceptably accurate at low frequency region (pitch acceleration term only provides high frequency contents).

By equating (2) and (5) on level ground, θ_(x)=θ_(xr), the following equation is obtained αX−βY=U  (11) where X=T _(roll)(z)(1−T _(ssc)(z)a _(y) Y=T _(roll)(z)(1−T _(ssc)(z)){dot over (ω)}_(x) U=T _(adi)(z)[ω_(x)+ω_(z)θ_(yr)]  (12)

The following is a summary of the computation algorithm.

In step 106, the filtered lateral acceleration X is determined in the following iterative equation X _(k) =d ₁₁ X _(k-1) +d ₁₂ X _(k-2) +d ₁₃ X _(k-3) +n ₁₁ a _(y) _(k) +₁₂ a _(y) _(k-1) +n ₁₃ a _(y) _(k-2) +n ₁₄ a _(y) _(k-3)   (13) where d_(1i),n_(1i), for i=1,2,3,4 are the involved first set of the filter coefficients.

In step 108, the filtered roll acceleration Y as is determined in the following iterative equation Y=d ₁₁ Y _(k-1) +d ₁₂ Y _(k-2) +d ₁₃ Y _(k-3) +n ₁₁{dot over (ω)}_(x) _(k) +n ₁₂{dot over (ω)}_(x) _(k-1) +n ₁₃{dot over (ω)}_(x) _(k-2) +n ₁₄{dot over (ω)}_(x) _(k-3)   (14)

In step 110 the filtered roll angle velocity, which is the sum of the roll rate and the product of the relative pitch and yaw rate, may be calculated in the following iterative algorithm U=d ₂₁ U _(k-1) +d ₂₂ U _(k-2) +d ₂₃ U _(k-3) +n ₂₁ R _(k) +n ₂₂ R _(k-1) +n ₂₃ R _(k-2) +n ₂₄ R _(k-3)  (15) where R_(k) is the total roll angle velocity at the kth time instant on level ground i.e., R _(k)=ω_(x) _(k) +ω_(z) _(k) θ_(yr) _(k)   (16) where d_(2i),n_(2i) for i=1,2,3,4 are the involved 2^(nd) set of the filter coefficients. If the approximation (10) is used, the total roll velocity on level ground can be approximated as in the following R _(k)≈ω_(x) _(k) +κω_(z) _(k) T _(pitch)(z)a _(x) _(k)   (17)

Using the calculated filtered value X_(k),Y_(k) and U_(k) at each sampling instant k, equation (11) can be used to gradient α, the roll acceleration coefficient β if the roll rate coefficient η is determined independent of the determination of α and β. A more complicated method is needed in order to obtain α, β and η simultaneously. One method to compute the unknown values α and β in equation (11) is the so-called least-square method.

Since equation (11) is true when the vehicle is driven on level ground and the vehicle does not have any wheel in the air (4 tires are contacting the road), a conditional least-square (CLS) method is used. Two CLS's may be used. The first CLS method conduct α and β parameter updating after a fixed number of conditional samples, while the second method updates α and β and resets a covariance matrix at each conditional sample.

Referring to FIG. 6A, the second conditional least square method is illustrated by box 250. Blocks 84, 86, 90 are the same as in FIG. 6. This embodiment, however, uses a relative pitch signal instead of the longitudinal acceleration signal. The relative pitch signal and yaw rate signal are multiplied in block 252. The output of block 252 is added to the roll rate in addition to block 95′ and filtered in block 96′. The filtered signals from blocks 86, 90, 96′ are provided to a conditional covariance block 254. The block 254 output is provided to roll gradient, roll acceleration coefficient least square updating block 256, which provides an output to block 100.

Since α and β are related to the inertia parameters of the vehicle body roll motion, only if the vehicle roll mode is fully excited α and β can be identified correctly through least-square methods. Hence, under the steady-state driving condition, the least-square identification is not conducted. Therefore, a further condition reflecting vehicle roll dynamic condition is needed. Considering the driver's steering input could induce roll dynamics, one of such conditions can be characterized using the steering wheel angle velocity, i.e., only if |{dot over (δ)}_(w)|≧δ_(wvmin) the CLS method will be conducted, δ_(w) is the measured steering wheel angle, δ_(wvmin) is a threshold (for example, 20 degree per second).

In step 114, the conditional sum of the products or cross products of the above filtered variables are determined over a significant large number N of the conditional samples. This is performed using the following iterative algorithm if the first CLS method is used

$\begin{matrix} {\begin{matrix} {{if}\mspace{14mu}\left( \mspace{14mu}{{{\overset{.}{\delta}}_{w}} \geq \delta_{{wvmi}\; n}} \right.} & {\&\mspace{14mu}{level}\mspace{14mu}{ground}} & \; \\ \; & {\&\mspace{14mu}{no}\mspace{14mu}{wheel}\mspace{14mu}{lifting}} & \; \\ \; & {{\&\mspace{14mu} s} \leq N} & ) \end{matrix}{\left\{ \mspace{14mu}{a_{11_{s + 1}} = {{a_{11_{s}} + {X_{k}^{2}\mspace{14mu} a_{21_{s + 1}}}} = {{a_{21_{s}} - {X_{k}Y_{k}\mspace{14mu} a_{22_{s + 1}}}} = {{a_{22_{s}} + {Y_{k}^{2}\mspace{14mu} b_{1_{s + 1}}}} = {{b_{1_{s}} + {X_{k}U_{k}\mspace{14mu} b_{2_{s + 1}}}} = {{b_{2_{s}} - {Y_{k}U_{k}\mspace{14mu} s}} = {s + 1}}}}}}} \right\}.}} & (18) \end{matrix}$

Notice that the time instant denoted by k is different from the CLS updating instant denoted by s. Only if the conditions are satisfied all the time, then s=k. N in (18) is the total number of conditional samples used for CLS, which could have value ranged from 1000 to 80,000.

Step 114 is performed for the conditional sample when the road is a level ground. The level ground can be identified if there is an indication that the vehicle is not on a significantly banked road. Therefore, checking road bank can be used for this purpose, for example, using the method disclosed in U.S. Pat. No. 6,718,248. The level ground can also be checked through a flatness index (as the one calculated in U.S. Pat. No. 6,718,248), or a road profile detection (see U.S. Pat. No. 6,718,248), or a rough comparison between the global roll angle and the nominal chassis roll angle.

In step 116, the nominal values for α₀ and β₀ of the roll gradient and roll acceleration coefficient are calculated. They can also be obtained by using test data from the vehicle with nominal loading condition and nominal configuration.

In step 118, the roll gradient α and the roll acceleration coefficient β are computed as in the following

$\begin{matrix} {{{{if}\mspace{14mu} s} = N}\left\{ \mspace{14mu}{{D = {\min\left( {{\max\left( {{a_{11_{N}} - \frac{a_{21_{N}}^{2}}{\max\left( {a_{22_{N}},a_{m\; i\; n}} \right)}},{- a_{m\; i\; n}}} \right)},a_{m\; i\; n}} \right)}};\mspace{14mu}{\alpha = {\min\left( {{\max\left( {{\left( {b_{1_{N}} - \frac{a_{21_{N}}b_{2_{N}}}{\max\left( {a_{22_{N}},a_{m\; i\; n}} \right)}} \right)D^{- 1}},\underset{\_}{\alpha}} \right)},\overset{\_}{\alpha}} \right)}};\mspace{14mu}{\beta = {\min\left( {{\max\left( {{\left( {\frac{a_{11_{N}}b_{2_{N}}}{\max\left( {a_{22_{N}},a_{m\; i\; n}} \right)} - \frac{a_{21_{N}}b_{1_{N}\;}}{\max\left( {a_{22_{N}},a_{m\; i\; n}} \right)}} \right)D^{- 1}},\underset{\_}{\beta}} \right)},\overset{\_}{\beta}} \right)}};\mspace{14mu}{s = 0};\mspace{14mu}{a_{11_{s}} = 0};{a_{21_{s}} = 0};{a_{22_{s}} = 0};\mspace{14mu}{b_{1_{s}} = 0};{b_{2_{s}} = 0};} \right\}} & (19) \end{matrix}$ where a_(min) is a small number (for example, 0.0001), which is used to remove division by zero in the implemented computation; α and α are the respective lower and upper bounds of the roll gradient, which can be calculated as α=α₀−Δα α=α₀+Δα  (20) and α₀ is the nominal value of the roll gradient (obtained through vehicle testing for the vehicle with nominal loading), Δα is the allowed variation of the roll gradient; β and β are the respective lower and upper bounds of the roll acceleration coefficient, which can be obtained as β=β₀−Δβ β=β₀+Δβ  (21) and β₀ is the nominal value of the roll acceleration coefficient (for vehicle with nominal loading), Δβ is the allowed variation of the roll acceleration coefficient. If the second CLS method is used, the roll gradient α and the roll acceleration coefficient β can be calculated in block 256 through the following iterative algorithm. First, a 2 by 2 matrix variable V_(s+1) at the (s+1)th conditional time instant is calculated from its past value V_(s) and the calculated filtered values of X_(k) and Y_(k) as in the following

$\begin{matrix} {\begin{matrix} {{if}\mspace{14mu}\left( \mspace{14mu}{{{\overset{.}{\delta}}_{w}} \geq \delta_{{wvmi}\; n}} \right.} & {\&\mspace{14mu}{level}\mspace{14mu}{ground}} & \; \\ \; & {\&\mspace{14mu}{no}\mspace{14mu}{wheel}\mspace{14mu}{lifting}} & \; \\ \; & {{\&\mspace{14mu} s} \leq N} & ) \end{matrix}\left\{ \text{}\mspace{14mu}{{{{if}\mspace{14mu}\left( {s = 0} \right)\mspace{14mu}\left\{ {}_{1}{{V_{0} = \begin{bmatrix} v_{0} & 0 \\ 0 & v_{0} \end{bmatrix}};} \right\}\mspace{14mu} V_{s + 1}} = {V_{s} - {g\;\frac{{V_{s}\begin{bmatrix} X_{k}^{2} & {{- X_{k}}Y_{k}} \\ {{- X_{k}}Y_{k}} & Y_{k}^{2} \end{bmatrix}}V_{s}}{1 + {{\gamma\left\lbrack {X_{k} - Y_{k}} \right\rbrack}{V_{s}\left\lbrack {X_{k} - Y_{k}} \right\rbrack}^{T}}}}}};} \right\}} & (22) \end{matrix}$ where g, γ and ν₀ are three positive numbers, and the time instant k denotes the regular time instant and the time instant s is the conditional time instant. Equation (22) is also called the covariance resetting in the normalized least square algorithm.

The roll gradient and the roll acceleration coefficient are calculated using the 2 by 2 matrix V_(s), and the calculated filtered values X_(k), Y_(k) and U_(k), as in the following

$\begin{matrix} {\begin{matrix} {{if}\mspace{14mu}\left( \mspace{14mu}{{{\overset{.}{\delta}}_{w}} \geq \delta_{{wv}\;\min}} \right.} & {\&\mspace{14mu}{level}\mspace{14mu}{ground}} & \; \\ \; & {\&\mspace{14mu}{no}\mspace{14mu}{wheel}\mspace{14mu}{lifting}} & \; \\ \; & {{\&\mspace{14mu} s} \leq N} & \left. \; \right) \end{matrix}\left\{ \mspace{14mu}{{\begin{bmatrix} \alpha_{s + 1} \\ \beta_{s + 1} \end{bmatrix} = {\begin{bmatrix} \alpha_{s} \\ \beta_{s} \end{bmatrix} - {g\frac{{S_{s}\left\lbrack {X_{k} - Y_{k}} \right\rbrack}^{T}U_{k}}{1 + {{\gamma\left\lbrack {X_{k} - Y_{k}} \right\rbrack}{S_{s}\left\lbrack {X_{k} - Y_{k}} \right\rbrack}^{T}}}}}};\mspace{14mu}{s = {s + 1}};} \right\}} & (23) \end{matrix}$

The calculated values are also limited to within their feasible sets as in the following α_(s+1) =sat(α_(s+1),α₀−Δα,α₀+Δα) β_(s+1) =sat(α_(s+1),β₀−Δβ,β₀+Δβ)  (24)

Notice that during the time instant where the conditions for (22) and (23) do not hold, the computations are frozen to the last values. That is if the following conditions are not true at a time instant k

$\begin{matrix} \begin{matrix} \left( \mspace{14mu}{{{\overset{.}{\delta}}_{w}} \geq \delta_{{wv}\;\min}} \right. & {\&\mspace{14mu}{level}\mspace{14mu}{ground}} & \; \\ \; & {\mspace{14mu}{\&\mspace{14mu}{no}\mspace{14mu}{wheel}\mspace{14mu}{lifting}}} & \; \\ \; & {\mspace{14mu}{{\&\mspace{14mu} s} \leq N}} & ) \end{matrix} & (25) \end{matrix}$ then computation in (2), (23) and (24) are not conducted, and the involved variables will keep the values for the last conditional sample s V _(k) =V _(s);  (26) α_(k)=α_(s); β_(k)=β_(s).

From the calculated roll gradient and the roll acceleration coefficient, the vehicle loading and its distance to the center of gravity of the vehicle body in the vertical direction can be determined in the following in step 122.

From FIG. 4, M_(s) is the vehicle body mass and the height of the center of gravity C₀ of the car body is denoted as h which is measured in the vertical direction with respect to the vehicle floor. A loading of mass ΔM is added to the vehicle. The distance between the center of gravity of this mass relative to the vehicle floor is denoted as H. The center of gravity C of the vehicle body is likely to change due to the added mass ΔM. The distance between C₀ and C is denoted as ΔH. Then ΔH must satisfy the following relationship

$\begin{matrix} {{\Delta\; H} = {\frac{\Delta\; M}{M_{s} + {\Delta\; M}}\left( {H - h} \right)}} & (27) \end{matrix}$ that is, the unknown ΔH can be readily obtained from the added mass ΔM, the height of the added mass H, the center of gravity height h of the vehicle body and the vehicle mass M_(s) prior to the addition of the mass ΔM.

The total roll inertial momentum of the vehicle body with respect to the final center of gravity C can be expressed as I _(xc) =I _(xc) ^(M) +I _(xc) ^(ΔM)  (28) where I _(xc) ^(M) =I _(xc0) +M _(s) ΔH ²  (29) I _(xc) ^(ΔM) =ΔM(H−h−ΔH)²

By inserting equation (27) into equation (29), equation (28) can be expressed as the following

$\begin{matrix} {I_{xc} = {I_{{xc}\; 0} + \frac{M_{s}\Delta\;{M\left( {H - h} \right)}^{2}}{M_{s} + {\Delta\; M}}}} & (30) \end{matrix}$

The nominal roll gradient and the roll acceleration coefficient is α₀ and β₀. Then

$\begin{matrix} \begin{matrix} {\alpha_{0} = \frac{M_{s}h}{K_{roll}}} \\ {\beta_{0} = \frac{I_{x}}{K_{roll}}} \end{matrix} & (31) \end{matrix}$ where K_(roll) denotes the roll stiffness due to suspension and anti-roll bar. Using those nominal values and the calculated values α and β, the loading mass and the loading distance satisfy the following

$\begin{matrix} {{{\Delta\;{MH}} = {{K_{roll}\left( {\alpha - \alpha_{0}} \right)} = A}}{\frac{\Delta\;{{MM}_{s}\left( {H - h} \right)}^{2}}{M_{s} + {\Delta\; M}} = {{K_{roll}\left( {\beta - \beta_{0}} \right)} = B}}} & (32) \end{matrix}$

From the relationship in equation (32), the following estimations may be determined

$\begin{matrix} \left\{ \begin{matrix} {H = {h + \frac{B}{2\; A} + \sqrt{{\frac{1}{4}\left( \frac{B}{A} \right)^{2}} + {h\left( \frac{B}{A} \right)} + \left( \frac{B}{M_{s}} \right)^{2}}}} \\ {{\Delta\; M} = \frac{\alpha}{h + \frac{B}{2\alpha} + \sqrt{{\frac{1}{4}\left( \frac{B}{A} \right)^{2}} + {h\left( \frac{B}{A} \right)} + \left( \frac{B}{M_{s}} \right)^{2}}}} \end{matrix} \right. & (33) \end{matrix}$

In step 124, the height H of the added mass (ΔM) from the vehicle floor based on the first equation of (33) is determined.

In step 126, the added mass Ml is determined based on the second equation of (33) by using the calculated height H of the added mass.

In step 128, a safety system such as a roll stability control system is controlled in response to the added mass and the height of the added mass. The safety system may also be controlled directly from the roll gradient and roll rate parameter, both of which may be adaptive. The safety system may also include the roll acceleration coefficient, which may also be adaptive as a factor in addition to the adaptive roll gradient and adaptive roll rate parameter. A vehicle roll angle compensation term may be determined such as that shown above. The vehicle roll angle compensation term may sensitize the roll angle computation for earlier deployments. The amount of roll compensation will depend on the particular vehicle and the characteristics of such things as the suspension and configuration of the vehicle. A roll stability control gain/threshold compensation term may also be generated. The threshold may be changed to allow earlier deployment if an added mass at a predetermined height is determined in steps 124 and 126. The amount of compensation will likely be determined experimentally based on the vehicle configuration.

Thus, as can be seen, the roll signal for control may be adjusted according to the vehicle loading and its height. On the other hand if the vehicle has a significant amount of loading, i.e., the calculated loading exceeds a threshold L_(max) ΔM≧L _(max)  (34) and at the same time the height of the mass exceeds another threshold H_(max) H≧H _(max)  (35) then in addition to the adjusted computation of roll signal for control (through the adjustment of the chassis roll angle using the new roll gradient and roll acceleration coefficient), all the gains for the feedback control used for commanding the actuators are set to a set of values G_(RLMAX) tuned for larger roof loadings, or adaptively adjusted in high values based on the amount of ΔM. In this way, when the vehicle has double wheel lift, the proper amount of actuation will be requested so the vehicle achieves robust roll stability control function during rollover event involving vehicles with large roof loadings.

If the vehicle has a significant loading, i.e., ΔM≧L_(max), but the loading is a floor loading, that is, the height of the loading is smaller than a lower threshold H_(min) H≦H _(min)  (36) all the gains for the feedback control used for commanding the actuators are set to a set of values G_(FLMAX), tuned for larger floor loadings.

If the vehicle has a significant loading, i.e., ΔM≧L_(max), but the loading height is between the roof and floor, i.e., H is between a lower threshold H_(min) (possibly zero) and the higher threshold H_(max) H _(min) ≦H≦H _(max)  (37) all the gains for the feedback control used for commanding the actuators are set to the following values adjusted based on the detected loading height as in the following

$\begin{matrix} {G_{H} = {G_{FLMAX} + {\frac{H - H_{\min}}{H_{\max} - H_{\min}}\left( {G_{RLMAX} - G_{FLMAX}} \right)}}} & (38) \end{matrix}$

If the vehicle has a roof loading which is below the maximum allowed roof loading L_(max) but above the lower bound of a valid roof loading L_(min), i.e., L _(min) ≦ΔM≦L _(max)  (39)

Assume all the nominal gains for feedback (for the vehicle with nominal loading) are denoted as G_(nom), then the control gains will be adjusted based on the detected roof loading as in the following

$\begin{matrix} {G_{\Delta\; M} = {G_{nom} + {\frac{{\Delta\; M} - L_{\min}}{L_{\max} - L_{\min}}\left( {G_{RLMAX} - G_{nom}} \right)}}} & (40) \end{matrix}$

The roll gradient itself can also be directly used to adjust control gains. If the vehicle has a significant increase of roll gradient, i.e., α≧α_(min), then all the gains for the feedback control used for commanding the actuators are set to the following values adjusted based on the detected roll gradient as in the following

$\begin{matrix} {G_{\alpha} = {G_{\alpha\;{MIN}} + {\frac{\alpha - \alpha_{\min}}{\alpha_{\max} - \alpha_{\min}}\left( {G_{\alpha\;{MAX}} - G_{\alpha\;{MAX}}} \right)}}} & (41) \end{matrix}$

where α_(min) is the roll gradient corresponding to a vehicle without roof loading and α_(max) is the roll gradient corresponding to the vehicle with maximum roof loading allowed.

Notice that the control gain adjustments other than the above listed linear interpolation methods are possible. Also, notice that the dead-bands and thresholds used in roll stability feedback control can also be similarly adjusted based on the loading height H and/or the loading ΔM, or the roll gradient, or the roll moment of inertia calculated as in (30).

Combining those control gain, dead-band and threshold adjustments with the quantitative roll angle (for example, the chassis roll) adjustment, a robust roll stability control function can be achieved.

Notice that based on the afore-mentioned computation of the added mass and its vertical from the vehicle floor, the roll moment of inertia of the vehicle body can be readily calculated as in (30). Alternatively, the roll moment of inertia I_(xx) can be determined directly by using the equation M _(s) a _(y) h _(cg) −K _(roll)θ_(xr) −D _(roll){dot over (θ)}_(xr) =I _(xx){dot over (ω)}_(x)  (42) where M_(s) is vehicle mass, a_(y), is the measured lateral acceleration, h_(cg) is the center of mass height above the roll center, K_(roll) is the roll stiffness of the suspension, θ_(xr) is the relative roll of the body with respect to the wheels (or ground), D_(roll) is the roll damping of the suspension, {dot over (θ)}_(xr) is the relative roll velocity and {dot over (ω)}_(x) the roll acceleration (obtained by differentiating the roll rate obtained from the roll rate sensor), and I_(xx) is the roll inertia, in kg-m².

When a_(y) and θ_(xr) are zero but {dot over (θ)}_(xr) and {dot over (ω)}_(x) are non-zero, then I _(xx) =−D _(roll){dot over (θ)}_(xr)/{dot over (ω)}  (43)

Or, when θ_(xr) and {dot over (θ)}_(xr) are zero and the mass, M_(s) is determined by another means, I _(xx) =M _(s) a _(y) h _(cg)/{dot over (ω)}_(x)  (44)

Alternatively, the roll response of the vehicle, obtained from the roll rate sensor, can be analyzed in the frequency domain to determine roll gradient and roll acceleration coefficient:

ω_(d)=the damped natural frequency of the roll motion of the vehicle (Hz)

ω_(n)=the undamped natural frequency of the roll motion of the vehicle (Hz)

=the damping ratio of the roll motion of the vehicle (unitless)

Taking the Laplace transform of

$\begin{matrix} {{\theta_{xr} + {\frac{D_{roll}}{K_{roll}}{\overset{.}{\theta}}_{xr}}} = {{\frac{M_{s}h_{cg}}{K_{roll}}a_{y}} - {\frac{I_{xx}}{K_{roll}}{\overset{.}{\omega}}_{x}}}} & (45) \end{matrix}$

and rearranging yields the following by considering {dot over (ω)}_(x)≈{umlaut over (θ)}_(xr) when the vehicle does not have large pitch motion or the vehicle does not have large yaw rate:

$\begin{matrix} \begin{matrix} {\frac{\Theta_{xr}(s)}{A_{y}(s)} = \frac{1}{{\frac{I_{xx}}{M_{s}h_{cg}}s^{2}} + {\frac{D_{roll}}{m_{s}h_{cgs}}s} + \frac{K_{roll}}{m_{s}h_{cgs}}}} \\ {= {\left( \frac{M_{s}h_{cg}}{K_{roll}} \right)\frac{1}{{\frac{I_{xx}}{K_{roll}}s^{2}} + {\frac{D_{roll}}{K_{roll}}s} + 1}}} \end{matrix} & (46) \end{matrix}$ where Θ_(xr)(s) and A_(y)(s) are the Laplace transform of the θ_(xr) and a_(y) respectively.

Using a digital fast Fourier transform to construct a frequency response function and determine the magnitude of

$\frac{\Theta_{xr}(s)}{A_{y}(s)}$

with respect to frequency, the roll gradient is equal to the DC gain (Equation 46 evaluated at s=0).

The peak of the magnitude of the frequency response function will occur at the damped natural frequency, ω_(d)=ω_(n)√{square root over (1−

²)}  (47) where ω_(n) is the undamped natural frequency and

is the damping ratio. (Alternatively, the damped natural frequency can be determined from “impulse-like” inputs into the suspension (hitting a bump in the road or a pot-hole), and determining the frequency of the roll oscillations). From (45), we have

$\begin{matrix} {\omega_{n} = {\left. \sqrt{\frac{K_{roll}}{I_{xx}}}\Rightarrow\frac{I_{xx}}{K_{roll}} \right. = {\frac{1}{\omega_{n}^{2}} = {- \beta}}}} & (48) \\ {\zeta = \frac{\omega_{n}D_{roll}}{2K_{roll}}} & (49) \end{matrix}$ Substituting Equation (49) into Equation (47) gives

$\begin{matrix} {\omega_{d} = {\omega_{n}\sqrt{\frac{\omega_{n}D_{roll}^{2}}{2K_{roll}}}}} & (50) \end{matrix}$

Equation (50) can be solved for ω_(n), which can in turn be used to evaluate Equation (48) and solve for the roll acceleration coefficient β.

Another way to determine roll inertia is to use an indirect empirical approach. First, a mass is determined, then the inertia may be correlated to mass. That is, a first estimate mass is determined via a method such as that described in Invention Disclosure (203-0482) filed as U.S. patent application Ser. No. 10/849,590, the disclosure of which is incorporated by reference herein, or another method. Then the inertia may be obtained per various methods explained below:

a. Assume inertia changes as mass changes from its base values, M_(s) _(_) _(base) to the current value, M_(s) _(_) _(current) as: I _(xx) _(_) _(total) =I _(xx) _(_) _(base) +h ₂(M _(s) _(_) _(current) −M _(s) _(_) _(base)   (51) where h₂ is assumed to be a constant. Based on testing for a vehicle loaded several ways, an average may then be used.

b. Correlate roll inertia, I_(xx), to vehicle mass; e.g., via testing using a vehicle inertial measurement facility to get roll inertia, the mass for various loading conditions and roof loads may be correlated. This is based on the premise that you can only load an SUV so many ways (e.g., mass above curb loading must be placed in vehicle seats and/or on the roof rack). Again, the vehicle mass may be determined via a method described in U.S. patent application Ser. No. 10/849,590 or another method, then use the roll inertia, I_(xx), corresponding to that mass obtained from testing.

The currently estimated mass or roll inertia estimate can be used several ways in the RSC algorithm, including:

Improved state estimation: Since mass and inertia values are used in the roll, pitch, and yaw angle state estimates, values accurately reflecting the actual vehicle will improve the accuracy, improve RSC control, and reduce the occurrence of false interventions.

Modified thresholds for activation: For example, a fully loaded vehicle can be less stable in roll. Hence, tighter thresholds on roll angle and roll rate, for example, can be used. This will increase the level of RSC control for the same level of roll angle and velocity for the loaded vehicle.

c. Increased gains for the RSC PID and/or Transition controllers.

d. Less stringent entry conditions for PID or Transition control intervention.

e. Any combination of the above.

Various vehicle dynamic and safety restraint characteristics may be adjusted upon the detection, determination, and evaluation of the loading or loading condition of a vehicle. The stated vehicle characteristics may be adjusted for any and all driving conditions and vehicle operating conditions. The vehicle characteristics may be adjusted for rollover and non-rollover events and situations. Several vehicle characteristics are described below in association with multiple embodiments of the present invention. The below described embodiments may be utilized alone, in cooperation with each other, or in any combination or order. Also, although the below embodiments are primarily described with respect to roof loading, they may be modified and utilized for other vehicle loading, such as interior vehicle loading, trunk loading, etc.

Throttle Adjustment

Upon the detection of a vehicle being lightly loaded, such as when there is little to no roof loading and/or when there is a single occupant and/or when there is little to no trunk loading, the throttle sensitivity and profile may be adjusted. The throttle sensitivity refers to the amount of throttle response per throttle adjustment. For example, the throttle response may be adjusted differently for the same amount of throttle displacement. The throttle profile refers to a throttle curve representing the relationship between throttle input and throttle response. In an example embodiment, the throttle adjustment is set such that the vehicle acceleration is not overly sensitive or is less sensitive in comparison with a vehicle that is heavily loaded. Upon the detection of the vehicle being heavily loaded, the throttle profile is adjusted to increase sensitivity. The adjustment is achieved such that the acceleration of the vehicle is more responsive in comparison to the vehicle without such loading.

Steering/Over-Steer/Under-Steer Characteristic Adjustment

When a vehicle has large roof loading, which is positioned close to the rear end of the vehicle, the vehicle tends to exhibit more over-steer than the same vehicle without such loading. Upon the detection of such loading, the operation of an electronic power steering (EPS) system may be adjusted to compensate for the experienced over-steer.

When a vehicle has large roof loading, which is positioned close to the front end of the vehicle, the vehicle tends to exhibit more under-steer than the same vehicle without such loading. Upon the detection of such loading, the operation of an electronic power steering (EPS) system may be adjusted to compensate for the under-steer.

When a vehicle has large roof loading, which is positioned close to the rear end of the vehicle, the vehicle tends to exhibit more over-steer than the same vehicle without such loading. Upon the detection of such loading, the operation of an active front steering system (AFS) may be adjusted to compensate for the over-steer.

When a vehicle has large roof loading, which is positioned close to the front end of the vehicle, the vehicle tends to exhibit more under-steer than the same vehicle without such loading. Upon the detection of such loading, the operation of an active front steering system (AFS) may be adjusted to compensate for the under-steer.

When a vehicle has large roof loading, which is positioned close to the rear end of the vehicle, the vehicle tends to exhibit more over-steer than the same vehicle without such loading. Upon the detection of such loading, the damping of the front suspension dampers of the vehicle is increased and the damping of the rear suspension dampers is decreased to compensate for the over-steer.

When a vehicle has large roof loading, which is positioned close to the front end of the vehicle, the vehicle exhibits more under-steer than the same vehicle without such loading. Upon the detection of such loading, the damping of the rear suspension dampers of the vehicle is increased and the damping of the front suspension dampers is decreased to compensate for the vehicle under-steer.

Vehicle Roll Motion Adjustment

When a vehicle has large roof loading, the vehicle tends to exhibit larger roll motions during dynamic maneuvers than the same vehicle without such roof loading. Upon the detection of such loading and when the vehicle is driven less aggressively, an EPS steering booster curve or associated profile may be adjusted to reduce the experienced driver head tossing in the roll direction. During aggressive maneuvers, the EPS booster curve is further reduced such that the roll motion of the vehicle is reduced for the same amount of large steering input, upon the detection of the roof loading.

Upon the detection of a roof loading of a vehicle, electronically controlled suspension dampers may be adjusted to adjust the vehicle roll damping characteristics. This occurs to reduce the increased vehicle roll motion due to the roof loading during vehicle cornering.

Vehicle Pitch Motion Adjustment

During braking, a vehicle with large roof loading experiences increased body-pitching velocity than the vehicle without such roof loading. Hence upon the detection of roof loading, an electronic brake distribution function with ABS may be adjusted such that the brake distribution increases more towards the rear axle. This reduces the vehicle body-pitching velocity, which further reduces driver head tossing in the pitch direction.

Upon the detection of roof loading, electronically controlled suspension dampers may be used to adjust the pitch damping of the vehicle such that the increased vehicle pitch motion is reduced during vehicle braking.

Safety Restraint Deployment Adjustment

When a roof of a vehicle is loaded, the acceleration, braking, steering, handling, response, and performance characteristics of the vehicle are different than without such loading. As such, safety restraints and the deployment thereof, including deployment rates, deployment profiles, and timing, may be adjusted to compensate for such characteristical changes. For example, if a vehicle has large, roof loading, its c.g. height increases. Hence when it is hit by another vehicle or its tire hits a tripped mechanism, the vertical distance between the impact spot and the vehicle's c.g. is likely to be large than the vehicle without large roof loading. Hence the vehicle will likely to rollover more than the vehicle without roof loading. For this reason, upon the detection of the large roof loading, the threshold for activating rollover curtain is adjusted.

Integrated Adjustment

Upon the detection of vehicle roof loading, the afore-mentioned individual adjustments through individual actuators may be arbitrated or prioritized among the various actuators so as to achieve a coordinated overall adjustment of the vehicle characteristics.

While particular embodiments of the invention have been shown and described, numerous variations and alternate embodiments will occur to those skilled in the art. Accordingly, it is intended that the invention be limited only in terms of the appended claims. 

What is claimed is:
 1. A method of controlling an automotive vehicle having an associated center of gravity, the method comprising the steps of: determining, in a controller, an added mass, onboard the vehicle in response to sensor signal information and relative to a known vehicle mass; determining, in a controller, a change in the associated center of gravity as caused by the added mass; controlling the automotive vehicle by adjusting at least one vehicle characteristic in response to the added mass and the change in the associated center of gravity.
 2. A method as in claim 1 wherein adjusting at least one vehicle characteristic comprises adjusting a vehicle braking characteristic.
 3. A method as in claim 1 wherein adjusting at least one vehicle characteristic comprises adjusting a vehicle steering characteristic.
 4. A method as in claim 1 wherein adjusting at least one vehicle characteristic comprises adjusting a vehicle throttle characteristic.
 5. A method as in claim 1 wherein adjusting at least one vehicle characteristic comprises adjusting a safety restraint deployment characteristic.
 6. A method as in claim 1 wherein adjusting at least one vehicle characteristic comprises adjusting a profile of the at least one vehicle characteristic.
 7. A method as in claim 1 wherein the at least one vehicle characteristic is adjusted in response to a comfort and ride handling disturbance.
 8. A method as in claim 1 wherein the at least one vehicle characteristic is adjusted in response to a roll condition.
 9. A method as in claim 1 further comprising: determining a loading condition in response to the added mass; and adjusting sensitivity of a throttle response in response to the loading condition.
 10. A method as in claim 1 further comprising: determining a loading condition in response to the added mass; and adjusting sensitivity of a steering response in response to the loading condition.
 11. A method as in claim 1 further comprising: determining a loading condition in response to the added mass; and adjusting sensitivity of a braking response in response to the loading condition.
 12. A method as in claim 1 further comprising: determining location of the added mass; and adjusting the at least one vehicle characteristic in response to the location.
 13. A method as in claim 1 wherein adjusting at least one vehicle characteristic comprises adjusting damping of a suspension.
 14. A method as in claim 1 wherein adjusting at least one vehicle characteristic comprises adjusting damping differently for each of a plurality of dampers.
 15. A method as in claim 1 wherein adjusting at least one vehicle characteristic comprises adjusting brake, distribution pressures differently for each of a plurality of brakes.
 16. A method as in claim 1 wherein adjusting at least one vehicle characteristic comprises adjusting a plurality of vehicle characteristics in a coordinated overall adjustment. 